Pulses of chaos synchronization in coupled map chains with delayed transmission 
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Pulses of synchronization in chaotic coupled map lattices are discussed in the context of transmis- 
sion of information. Synchronization and desynchronization propagate along the chain with different 
velocities which are calculated analytically from the spectrum of convective Lyapunov exponents. 
Since the front of synchronization travels slower than the front of desynchronization, the maximal 
possible chain length for which information can be transmitted by modulating the first unit of the 
chain is bounded. 
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I. INTRODUCTION 

Chaos synchronization is a counter-intuitive phe- 
nomenon which has been extensively investigated since 
its original discovery [H, [!]• In particular, its potential 
of being applied for novel secure communication devices 
has attracted a lot of research on the foundations of chaos 
synchronization . Motivated by experiments on syn- 
chronized chaotic lasers, networks of chaotic units with 
time-delayed couplings and feedback are recently being 
studied [y, 0] ■ For chains of chaotic units with delayed 
couplings other counter-intuitive phenomena have been 
found such as anticipated or sublattice synchronization 
[1, Q . Chains of chaotic units are also discussed in the 
context of convective instabilities. Instabilities of such 
systems have successfully been described by spectra of 
comoving Lyapunov exponents [Tol - [l5| . 

There are several possibilities to transmit a message 
between two synchronized chaotic units A and B [IGI]. 
One possibility is chaos modulation: the message modu- 
lates the dynamics of unit A. This affects the dynamics 
of unit B and it is possible to recover the message with a 
certain bit error rate [U [3 • This application is the mo- 
tivation for the following question which we investigated 
in this report: is communication based on chaos modu- 
lation possible in a long chain of unidirectional coupled 
chaotic units? 

The dynamics of the first unit is changed according to a 
given bit sequence. These bits are recovered by observing 
the synchronization of the two last units of the chain. 
Note that we keep the system simple: the bits are not 
recovered for each unit but only with the two last units 
of the chain indicated in Fig. [T] The intermediate units 
act as passive relay stations which may be necessary to 
transmit signals over long distances. Note that the units 
can transmit the signals with arbitrary delay times. In 
order to detect perturbations of the first unit at the end of 
the chain we introduce a self-feedback of the first sending 
unit with a delay time which is identical to the coupling 
delay of the last two units. We analyze the transmission 
of information using the method of comoving Lyapunov 




FIG. 1. Scheme of the examined chain: each unit is driven by 
its predecessor. 



exponents. 

In fact already for single units with delayed feedback 
one finds, in a spatiotemporal representation, convective 
instabilities which move in a cone with a spectrum of ve- 
locities [l^. A chain of oscillators, nonlinear partial dif- 
ferential equations close to instabilities, and lasers with 
delayed feedback have been analyzed with the spectrum 
of comoving Lyapunov exponents . Chains of chaotic 
units which have been discussed in the context of turbu- 
lence show convective instabilities as well , even if 
the single units have a delayed self-feedback which gen- 
erates anticipated chaos jl5l |. 

In this report we extend this analysis to time-delayed 
couplings in the context of secure communication. 



II. SYSTEM DYNAMICS 

In this report we investigate the previous general ques- 
tion for a simple model: a chain of coupled maps with 
time-delayed transmission. The first unit has a self- 
feedback which is necessary to obtain zero lag synchro- 
nization for the last two units of the chain. The system 
is defined by the following equations: 

x° = (l-e) f{xU)+e f{xtr,) 

x- = {l-e) f{xU)+e f{x-zl) (1) 

is the variable of the nth unit at the discrete time 
step t. f{x) is a map which yields chaotic iterations. The 
transmission time from unit n — 1 to unit n is denoted by 
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Tn and e is the coupling strength. The first unit n = 
has self- feedback with delay tq. It is easy to see that the 
system has the following synchronized solution: 



A„ = ^ Tfc - n ■ To (2) 

fc=i 

If all delays are identical to the feedback delay of the 
first unit, r„ = tq, one finds complete zero lag synchro- 
nization, x" = Xf. For smaller values of r„ one finds 
anticipated chaos whereas for larger ones the units lag 
behind the first one. When the transmission delay from 
unit iV — 1 to unit N is identical to the feedback tq of 
the first unit, the last two units are synchronized with- 
out time shift, x^~^ = xf . Hence, the synchronization 
of the chain can be measured from the difference of the 
variables of the last two units. From now on identical 
delay times t„ — tq are used, since different r„ lead to 
trivial time shifts. The stability of the chaotic synchro- 
nized solution x" — x^ can be calculated by considering 
small perturbations of this trajectory and linearizing Eq. 
([1]). The evolution of deviations between two consecutive 
units, J" = — a;"~^, n = 1, . . . , TV are determined by 
the following set of linear equations: 

5l = {l-e) f\xU)5U 
5- = {l-e) nxU)5U + e fixU^Klr, (3) 

If only the unit n is perturbed, 5^ — for k < n, then 
the deviation follows the equation 

= {l-e) nxU)5U (4) 

This perturbation is stable if 

(ln|/'(x?)|)<-ln(l-e) (5) 

Note that the average (. . .) is taken over the first unit 
which, due to the feedback delay tq, is hyperchaotic; i.e., 
it has a spectrum of positive Lyapunov exponents. Ac- 
cording to Eq. ([5]), the synchronized solution ([2]) is stable 
if the coupling strength e is large enough. 

III. CHAOS MODULATION 

Now we apply chaos modulation to a concrete system: 
we take f{x) — rx{l — x) with r = 4 for the bulk of 
the chain and a variable rg for the first unit depending 
on the message to transfer. If rp = r = 4 then the 
chain will successively relax into the synchronized state 
starting from the head. If one chooses tq = r' ^ r, 
for example, r' = 3.95, the system leaves the common 
trajectory, again beginning at the head. This creates 
pulses of (de) synchronized phases traveling throughout 
the chain. We will call rg = r the tuned and ^ r 
the untuned state. Figure [5] shows numerical results of 
a simulation for e = 0.8, tq = 10: the chain starts in a 
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FIG. 2. Absolute value of the difference 5" between neigh- 
boring units: for t < 0, 100 < t the first unit is un- 
tuned (ro = r' = 3.95). For < t < 100 the chain is 
tuned (ro = r = A), thus, creating a synchronized pulse. 

£ = 0.8, TO = 10. 

desynchronized state and is tuned at t = 0. After 100 
time steps it is untuned again. The absolute difference 
(5" I = \xt — is shown in color code; i.e., the dark 

areas mark synchronization. 

The velocities of the signal fronts are defined as 

Anumber of units 

^ = 

Atime 

In numerical simulations two units are considered syn- 
chronized if |(5"| < O with a small threshold 8 « 10"^". 
We find that for <C 1 the results for the velocities of 
the head and tail of the pulse do not depend on the actual 
choice of the threshold. In the simulation shown in Fig. 
[5] one finds that synchronization travels with a speed of 
approx Vs = 0.088 whereas desynchronization propagates 
with Vd — 0.100 which is the maximum possible velocity 
defined by the coupling delay tq = 10. Thus, desynchro- 
nization finally overtakes synchronization and the pulse 
cannot be detected in the rest of the chain. 

Figure[3]shows Vg as a function of the coupling strength 
e. Near the critical coupling one observes a slowing down. 
One would expect a monotonous decrease in the speed 
with decreasing coupling strength. However, there are 
two small peaks in the data. For tq > 1 the value dis- 
tribution of the first unit and thus the local Lyapunov 
exponent Aq = (ln|/'|) averaged over this distribution 
are functions of e. Near the peaks Aq becomes rather 
small making the map "less chaotic" or even non-chaotic 
and thus accelerating synchronization. 

Figure |4] shows the additional delay per unit 7 with 
z; = 1/(to +7) as a function of tq. For large tq there is a 
saturation effect. 



IV. CONVECTIVE PERTURBATIONS 

For a chain of connected relay stations one might ex- 
pect the coupling delay between the units to be on a 
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FIG. 3. Propagation velocity of synchronization in units of 
the maximum velocity 1/to for different coupling parameters 
e at To = 10. Here the critical coupling is Ec ~ 0.35. 
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FIG. 4. The additional delay 7 — 1/vs — ro as a function of 
the coupling delay ro with e — 0.6: for large ro there is a 
saturation effect. 

much larger time scale than the one of the internal dy- 
namics. In our model this would imply tq 3> 1. However 
the numerical results presented by now indicate that the 
value of tq has a quantitative effect on the velocities but 
does not change the behavior of the chain in principle 
(see Fig. d]). Hence, for analytical treatment we shall fo- 
cus on To = 1. We start with a completely synchronized 
chain, = x^, and consider small perturbations, Eq. 
Then the discrete Green's function of the system is 

5*"= (^) ii~eY-e"f[f(xl,) (6) 

The evolution of an arbitrary perturbation Sq of the syn- 
chronized solution can be obtained by 

k 

The Green's function can be approximated by 

\g^\^eMmt), v = n/t (7) 

A(t;) = (l-i;)ln(^i-|^+t-ln(^)+Ao (8) 
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FIG. 5. The convective Lyapunov exponent A(ii), Eq. ([SI, 
for Ac = 0.3 and e — 0.6. 

where Aq = {\n\f'(xf)\) is the Lyapunov exponent of 
the applied map and A(u) is the convective Lyapunov 
exponent as mentioned in [lT| - [T5| . Outside the interval 
w G [0 : 1] A(w) has to be set —00 due to causality. Figure 
[5] shows the convective Lyapunov exponent A(w) of Eq. 
dH]). A{v 0) describes the relaxation of a local pertur- 
bation, for negative values the chain finally relaxes to the 
synchronized trajectory, in agreement with Eq. The 
fronts of desyncronization and resynchronization are at 
the border of stability, A{v) = 0. Hence, the first zero 
at Vs determines the velocity of synchronization and the 
second one yields the velocity of desynchronization. If 
ln(e) > — Ao one finds the maximal possible desynchro- 
nization speed W(j = 1. 

Hence, a local perturbation of the synchronized tra- 
jectory generates a pulse of deviations with a front ve- 
locity and a tail velocity Vg. The maximum of the 
pulse propagates with Vm = s and the integrated pulse 
increases exponentially with the exponent Aq. 

Figure [6] compares this theory with results obtained 
from simulations. Obviously the prediction of A(w) is 
correct even if the linear equations do not describe the 
dynamics in the noisy regions of the chain preceding the 
final synchronization process. Note that for tq = I one 
finds Ao = In 2; hence, the critical coupling is Sc — 1/2. 
At this critical point the velocity of synchronization goes 
to zero; again one observes a critical slowing down for 
synchronization. 



V. INFORMATION TRANSFER 

Pulses of synchronization and desynchronization can 
be used for communication: sender and receiver define a 
certain bit length L. Then the sender switches tq between 
r and r' in time intervals L according to the bit sequence 
(let us define tq = r = true, vq = r' = false). Regions 
of corresponding state are created at the head and start 
to propagate throughout the chain. The receiver decodes 
the message by observing the synchronization state of the 
last two units. 
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FIG. 6. Vs and Vd obtained by simulations are compared to 
the predictions made by h.{vi) — for tq = 1. 



It is in any case Vg < Vd so the true-bits will always suf- 
fer from shrinking and a single true-bit will finally vanish 
at a chain position approximately determined by 



L ■ 



Vd Vs 
Vd - Vs 



(9) 



than isolated ones such as 01010. 



VI. CONCLUSION 

Pulses of synchronization in chaotic chains can be used 
for communication. The information has not to be re- 
stored at every unit but is recovered by observing the 
synchronization status of the last two units. We have 
analyzed the propagation of pulses using the method of 
convective Lyapunov exponents JlO-15, 19). We find that 
the front of a synchronization pulse travels slower than 
its tail. This limits the maximal chain length for which 
the pulse can be detected and, thus, has a large influence 
on the bit error rate. 

We have investigated the propagation of binary infor- 
mation by detuning the first chaotic unit. For a chain of 
chaotic lasers this may be accomplished by detuning the 
pump current of the first laser. It would be interesting 
to investigate the conditions for which analog signals can 
be transmitted. From our results we expect two major 
criteria for reliable transmission: The signal must have 
a spectrum of low frequencies compared to the inverse 
minimal bit length ^ and the linear approximation of 
the time evolution must hold; see Eq. ([3]). 



However if a true-bit is followed by a second one, the 
second bit can travel a longer distance before vanishing. 
Obviously bit sequences with consecutive ones such as 
01110 can be transmitted with much lower bit error rates 
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